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Quantum nonlocality and contextuality are two phenomena stemming from nonclassical correla¬ 
tions. Whereas the former requires entanglement that is consumed in the measurement process the 
latter can occur for any state if one chooses a proper set of measurements. Despite this stark differ¬ 
ences experimental tests of both phenomena were similar so far. For each run of the experiment one 
had to use a different copy of a physical system prepared according to the same procedure, or the 
system had to be brought to its initial state. Here we show that this is not necessary and that the 
state-independent contextuality can be manifested in a scenario in which each measurement round 
is done on an output state from the previous round. 


Introduction. In a standard Bell scenario a pair of ob¬ 
servers share a bipartite system on which they perform 
local measurements mm- The results of these measure¬ 
ments may not be explainable by local realistic theories, 
however to observe it one needs a quantum system pre¬ 
pared in an entangled state. The entanglement contained 
in this state is consumed during the measurement and 
the resulting post-measurement state is local and useless 
for further Bell tests. A similar effect, although more 
subtle, takes place in the state-dependent contextuality 
scenarios. For example, in the Klyachko-Can-Binicioglu- 
Shumovski (KCBS) test |5] an initial state of the system 
can exhibit contextuality with respect to a specific set of 
measurements, but all the post-measurement states are 
noncontextual if tested in the same KCBS test. 

It is therefore natural to think of states exhibiting 
nonlocality or contextuality as some resourceful states, 
whereas the remaining states can be considered as re¬ 
sourceless. In this sense, the resourceful states can pass 
the test at the cost of becoming resourceless. Since every 
such test requires a sufficient amount of data to statis¬ 
tically determine its outcome, more than one measure¬ 
ment has to be performed. This requires an ensemble of 
resourceful states from which one draws a system in each 
measurement round, or a resetting procedure in which 
one brings back resourcefulness to the post-measurement 
state. 

The above interpretation makes the state-independent 
contextuality a different phenomenon. Every quantum 
state of a more than two-level system can exhibit contex¬ 
tuality if one prepares a special set of measurements |4j- 
Ej. Due to this fact one cannot divide the set of all states 
into resourceful and resourceless since there is no resource 
consumption. Therefore, one is inclined to ask: How to 
reuse post-measurement states in some state-independent 
contextuality scenario ? 
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There are two additional motivations behind this ques¬ 
tion. First of all, if there were an efficient method of 
state-recycling it would radically simplify any experi¬ 
mental implementations of contextuality tests in which 
measurements are non-destructive (e.g. trapped ion ex¬ 
periments 0)- Moreover, from the fundamental point 
of view it is still unclear what kind of resource contex¬ 
tuality is and how to quantify it with respect to some 
meaningful tasks jlOj . Showing that post-measurement 
states can be efficiently reused in state-independent con¬ 
textuality tests would suggest that in order to look for 
resourcefulness of contextuality one should not associate 
it to nonlocality and entanglement. A search for a new 
fundamental meaning of contextuality is in order. 

In this work we propose a state-recycling method to 
investigate contextuality. We show that in the Peres- 
Mermin state-independent contextuality scenario EH3, 
which can be also expressed in a form of an inequality 0, 
the post-measurement states can be reused in the next 
measurement round. Moreover, because the set of post¬ 
measurement states is finite, if sufficiently many measure¬ 
ment rounds are performed each state will be measured 
in every measurement context. As a result, the measure¬ 
ment process on recycled states can be described as a 
Markov chain from which one obtains necessary statis¬ 
tics to evaluate the result of the test on all the states. 
Alternatively, it can be viewed as a test on an effective 
state that corresponds to the stationary distribution of 
the process. Eventually, we consider imperfect measure¬ 
ments and show that even in the presence of imperfect 
measurement settings our state-recycling model is still 
capable of detecting contextuality. 

State-independent contextuality scenario. Contextual¬ 
ity is a general phenomenon that can be formulated out¬ 
side of quantum theory and as such it can be studied 
within the operational framework of black boxes nnna. 
This framework uses only the concepts of preparation, 
transformation and measurement. However, we want 
to study contextuality in realistic systems and the only 
known contextual systems are quantum ones. We there- 
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fore study the problem of state-recycling within the quan¬ 
tum formalism. 

We consider the scenario commonly known as the 
Peres-Mermin square EK2- It consists of nine dichotomic 
±1 observables on a four-level system for which one can 
distinguish two two-level degrees of freedom 

A± = a x (g) 1 A 2 = 1 0 <7 y A 3 = a x G> a y 

A 4 = 1 ig) a x A 5 = a y ® 1 A e = a y (g) a x (1) 

At = cr x ® crx A s = a y <g> <J V A$ = cr z ® a z 

where cq are Pauli matrices. Each row and column con¬ 

sists of three jointly measurable observables. The as¬ 
sumption of noncontextuality states that a measurement 
outcome of each observable does not depend on the mea¬ 
surement context, i.e., whether it is measured with the 
observables from the same row, or with the ones from the 
same column. 

However, the above assumption is invalid because 
quantum observables are contextual. Triples of the ob¬ 
servables in each row and column are jointly measurable 
because of a pairwise commutation. Note that the prod¬ 
uct R 3 of the three operators in the last row yields — 1 
whereas the products (Ri, i = 1,2) of three operators in 
any other row or column (C), i = 1,2,3) give +1. This 
gives riie{i 2 3} -RjCi = —1- However, in order to calcu¬ 
late this product each operator is used twice. Therefore, 
noncontextuality implies that no matter what value is 
assigned to each observable this product should be equal 
to -f*l- 

The above paradox can be expressed in terms of an 
inequality |8] 

(A 1 A 2 A 3 ) + (A^A^Aq) — {AjAgAg) 

+ {A\A 4 At) + {A 2 AgAg) + {AgAgAg) < 4 (2) 

The upper bound for this inequality comes from a simple 
optimization procedure over all possible values Ai = ±1 
that takes into account the noncontextuality assumption. 
It is straightforward to show that for the quantum spin 
operators we will find the left-hand side of ([ 2 ]) equal to 
6 because the products of all triples are 1 except for 
AyAgAg = — 1. 

Measurements. There are six possible measurement 
contexts in the Peres-Mermin scenario and in principle 
the measurement of each context can give one of eight 
possible outcomes: + + +, + H—, etc. However, if the 
measurements are perfectly implemented, there can be at 
most four different outcomes. This fact lies at the very 
root of quantum state-independent contextuality. In this 
case, a measurement of some triple of observables is a 
projection of the quantum state of the system onto one 
of four eigenstates of the triple, irrespective of whether 
the triple of measurements were done simultaneously or 
sequentially. 

Let us define the triple-eigenstates as b lt ), where 
j = 1, ...,6 denotes the triple and i = 1, ...,4 denotes the 


(a) (b) 



FIG. 1. Schematic representation of the standard test of con¬ 
textuality in which each measurement round is performed on 
a different system drawn from an ensemble (a) and of the 
state-recycling scenario in which each measurement round is 
performed on the same system (b). 


basis state. There are 24 different states, 16 local (cor¬ 
responding to tensor products of eigenstates of a x and 
a y - triples j = 1,... ,4) and 8 nonlocal (4 Bell states 
corresponding to the last row j = 5 and 4 rotated Bell 
states corresponding to the last column j = 6). 

State-recycling scheme. In the standard contextuality 
scenario one performs each measurement round on a dif¬ 
ferent system drawn from an ensemble. It is assumed 
that the ensemble is described by the state p. For each 
measurement round an experimenter randomly chooses 
one of six measurement contexts and after many rounds 
one obtains enough data to evaluate the inequality © 
see Fig. [I] (a). 

Here we propose a different approach in which a single 
copy of a system, initially prepared in some state po, is 
measured in a randomly chosen basis j and as a result 
ends up in some post-measurement state | bj.) called a 
triple state. Without loss of generality we can assume 
that po is one of the 24 triple states. Then, the experi¬ 
menter randomly chooses which triple j' to measure next 
and the same system is plugged again to the measuring 
device. The system ends up in the triple state | by,). This 
procedure is repeated for N rounds - see Fig. 0(b)- 

In each round the probability for a state | by) ending 
up as | by) is given by 

Tud', )> ( 3 ) 

where Mj. = The random choice of the mea¬ 

surement triple gives us the factor of 1/6. Note that 
the quantity Ty y is symmetric, i.e., Tj. y = Ty y and 
that its evaluation requires no knowledge of any other 
earlier state. Therefore, the many-round measurement is 
a Markov chain defined on a 24-state space. 

The state of the system after t measurement rounds 
can be represented as a 24-dimensional probability vector 
p(t) describing a probability distribution over all triple 
states. Note, that although the triple states are not all 
mutually orthogonal, we can distinguish the outcomes of 
measurements because we store information about which 
triple is measured and the four eigenstates of each triple 
are mutually orthogonal. 





3 


The transition matrix T of a process is given by Eq. 
([3]) - see Appendix for an explicit form. One step of this 
process is given by p(t + 1) = Tp(t). It is easy to verify 
that after many rounds the state of the system converges 
to the stationary distribution n = Tir which is a uniform 
distribution over all 24 triple states. 

Interpretation. We consider two alternative interpreta¬ 
tions of the contextuality test based on the above Markov 
chain. In the first interpretation the state-recycling 
scheme leads to tests of contextuality on all 24 triple 
states. Note that if the process starts in a state b lt ), then 
after sufficiently many rounds it comes back to the same 
state. A crucial observation is that physically the state 
\bj x ) that is obtained via recurrence has exactly the same 
properties as the initial state. This leads to a conclusion 
that from the experimental point of view one cannot ob¬ 
serve a difference between state preparation in a scenario 
in which subsequent states are drawn from an ensemble 
and a scenario in which subsequent states appear as a 
result of the recurrence in a state-recycling scheme. 

After the recurrence the next measurement is chosen 
randomly, therefore for a large number of rounds one can 
obtain enough data to evaluate the inequality @ f or I&*) 
as an input state. However, since we already know that 
the system state converges to the stationary distribution 
7 r, which is evenly distributed over all possible states, one 
can evaluate the inequality ([ 2 ]) for any triple state, given 
that sufficiently many measurements were performed. 

The expected number of steps to achieve the stationary 
distribution 7r with an accuracy e, i.e., the mixing time, 
is bounded from above by 

, . 3, /24\ 

tmix(e) < -log —J . (4) 

See Appendix for details. For e of the order of 10 -3 , 10 —5 , 
and 10 -10 the mixing times are less than 16, 23, and 
40, respectively. If we further assume that one gathers 
data once the system reached the distribution 7r then 
each triple state can be obtained with equal probability 
A. One can show that an average number of sequential 
measurements one needs to perform on the stationary 
state to detect all triple states is 91 - see Appendix. 

Alternative interpretation. Instead of considering con¬ 
textuality tests on all 24 triple states, one can consider a 
single test of contextuality on a stationary state n. Be¬ 
cause of indistinguishability between a stationary state 
and an ensemble represented by the quantum state 

= \ b jt)( b 3i\ = \’ ( 5 ) 

ij 

measurement data collected after | log (A j steps can be 
considered as data measured on the state p n . 

There is however one disadvantage in the above reason¬ 
ing. One has to be aware of the fact that after the mea¬ 
surement the state is bj i ), therefore the next measure¬ 
ment is not performed on p However, the knowledge of 


the triple state that is plugged into the next measuring 
device is conditioned on the knowledge of the previous 
measurement outcome. In the end, the averaging pro¬ 
cedure that is applied to evaluate <@> effectively reduces 
the post-measurement state to p n . 

Note, that in a sense this problem is also present in the 
standard scenarios in which an ensemble is represented 
by a mixed state The only difference is 

that in our case the knowledge of an exact preparation 
is known to the experimenter, whereas in the ensemble 
case this knowledge is encoded on some other inaccessible 
system. Therefore, the above interpretation relies on the 
assumption that this difference is irrelevant. In addition, 
in the standard ensemble scenario one can try to solve 
this problem by preparing a larger ensemble of entangled 
systems JA In this case one can use Bell- 

type arguments that the knowledge of the sub-ensemble 
preparation simply does not exist prior to a measurement 
on an auxiliary system. 

Errors. While performing a sequence of measurements, 
we cannot ignore errors which are created by fluctuations 
of control parameters m These errors will lead to a 
detection of more than four outcomes and hence the in¬ 
equality <[2j) is no longer maximally violated [14j . The 
errors occur because of a lack of control and in particu¬ 
lar if measurement device settings fluctuate one cannot 
determine the post-measurement state with a perfect ac¬ 
curacy. 

In our error model we assume that errors come from an 
imperfect alignment of measuring devices. Each triple of 
measuring devices is perfectly aligned with a probability 
p and completely out of control with a probability 1 — p. 
Therefore, the four outcomes corresponding to the triple 
states are detected with the probability p and the remain¬ 
ing four outcomes, to which we refer as error states , are 
detected with the probability 1 — p. 

From the point of view of the Markov chain describ¬ 
ing the imperfect scenario, the probability space is now 
48-climensional. To calculate a new transition matrix we 
assume that each error state is represented by the maxi¬ 
mally mixed state . Note that even if in some measure¬ 
ment round an error state is detected, in the next round 
one can still detect a triple state with the probability p. 
Therefore, the state-recycling scheme is self-correcting. 

The new transition matrix can be schematically repre¬ 
sented in the form 


T\ 

t 2 

t 3 

T 4 


where T\ = pT is the previous matrix Ml) multiplied by 
p. The other parts of the above matrix Ml) refer to tran¬ 
sitions to, from and between the error states. They read 
T 2 = 24 11-24x24 and T 3 = T 4 = ^1 2 4x 24, where t nxn 
is the n x n matrix with all elements equal to one. In 
brief, the matrix T 2 describes the transitions from the 
error states to the triple states. The opposite process 
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is encoded by T 3 . Finally, T 4 describes transitions be¬ 
tween the error states. Note that T err is stochastic as its 
columns sum up to 1. 

The stationary distribution of T err corresponds to 
its eigenvector with the eigenvalue 1 and is ir = 
(m> • • •) 2 ii ~^ai • • • > 33 ^), where the first 24 entries cor¬ 
respond to the triple states and the last 24 to the error 
states. It means that after sufficiently many steps all 
triple states are equally probable to be detected. This is 
also true for the error states. Moreover, the probability 
of the triple state detection is the same as the probability 
p of the perfect setup alignment. 

In the presence of noise the violation of the inequality 
([ 2 ]) is not maximal. If we assume that the data is collected 
after the stationary state is obtained, the influence of 
noise can be easily calculated. Since the triple states 
appear with the probability p and the error states with 
the probability 1 —p, one has (AiAjAk ) = 2 p—l (1 — 2 p 
in the case of {A 7 A$Ag)). As a result, the inequality ([2| 
reads 

12 p-6< 4. (7) 

It leads to the conclusion that the accuracy of the mea¬ 
surement setup has to be greater than p > 5/6 « 0.83 to 
see the violation. 

Discussion. Despite the advantage of using a single 
copy of a system, the state-recycling scheme can be vul¬ 
nerable to a memory loophole. Apart from the loophole 
already discussed in Ref. m, it might be possible that 
the system stores information about the measurement 
triples measured in the past allowing for some contex¬ 
tual hidden variable model. This model can take advan¬ 
tage of this data to mimic the contextual behaviour. One 
can reasonably assume that in scenarios in which every 
measurement round is performed on a different system 
drawn from an ensemble such a possibility is highly un¬ 
likely since the stored information would have to be com¬ 
municated from one system to the other. However, in the 
state-recycling scheme every measurement is performed 
on a single system, therefore the contextual hidden vari¬ 
able model would not need additional communication. 

Although we do not construct any such model, we ar¬ 
gue that its existence would require additional resources 
m- In the state-recycling scheme one tests a quantum 
system on which only two bits of classical information can 
be efficiently stored. On the other hand, to store infor¬ 
mation about which triple was measured one requires an 
additional log 2 6 bits of an auxiliary memory per triple. 
Therefore, in order to simulate contextuality on a classi¬ 
cal system one would require a system capable of storing 
more data than the original quantum system. 

One has to be aware that the problem of hidden vari¬ 
ables in contextuality tests is slightly different from the 
one in Bell tests. In the case of Bell tests there exists a 
strong physical argument against the possibility of local 
realistic description (assuming the free will of observers). 


In order to reproduce the results of the Bell test on classi¬ 
cal systems one would require superluminal communica¬ 
tion, which is forbidden by the special relativity theory. 
On the other hand, in the case of contextuality there are 
no strong physical arguments against contextual hidden 
variables. The arguments against them are of practi¬ 
cal nature, i.e., a simulation of contextuality on classical 
systems would be inefficient from the point of view of 
resources m- Following Occam - it is more efficient to 
describe contextuality as a nonclassical effect. 

Finally, we would like to stress that our motivation 
is not to propose a test of contextuality that is more 
robust against loopholes. We propose a test whose im¬ 
plementation is fundamentally different and requires less 
experimental preparation. 
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APPENDIX 
Transition matrix 

Here we give an explicit form of the transition matrix 
governing the state-recycling scheme 
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It is symmetric and one can verify that the eigenvector 
corresponding to the eigenvalue 1, the stationary state, 
is of the form 

n = ( 0 ) 

The remaining eigenvalues are ^ (9 eigenvalues) and 0 
(14 eigenvalues). 

Mixing time 

Next, we calculate the expected number of steps to 
achieve the stationary distribution ir (we follow Ref. 
m)- First, we need to define a distance between two 
probability distributions p and v on a finite set id 

Wli-vW =max\n(A)-v(A)\. (10) 

In our case = {|6j.)}. Next, we define a quantity which 
tells us how far from stationary probability we are in the 
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t -tli step 

d(t) = max ||p(f) g — 7r g ||, (11) 

where q denotes the g-th coordinate of a probability vec¬ 
tor. Then the mixing time, i. e., time t needed for p(t) 
to be approximated by 7r, is given by 

imix(e) = min {t : d(t) < e} , (12) 

where e denotes the accuracy of the approximation. 

For the above Markov chain the following inequality 
holds 


where 7r m j n = min g 7 T q and A* < 1 is the second greatest 
eigenvalue of T. Since A* = | we get 

imix(e) < ^ log . (14) 
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Expected number of measurements to detect all 
triple states 

We use the solution of coupon collector’s problem to 
calculate the average number of times we need to keep re¬ 
peating the measurement procedure to detect all n = 24 
possible triple states. Originally, the coupon collector’s 
problem concerns how many times on average one needs 
to draw a coupon with replacement from some larger en¬ 
semble to collect all different coupons. There are many 
different generalisations of this problem, e.g. to the sce¬ 
narios of non-uniformly distributed probabilities of find¬ 
ing each coupon or collectors not necessarily looking for 
all possible coupons m- 

Here we consider a general case of the state-recycling 
protocol in which errors can occur. We can formulate our 
problem of finding all possible triple states of the system 
in the state-recycling scenario as collecting coupons with 
almost-uniform probability distribution, meaning that an 


extra coupon, which does not belong to the collection 
(error state), might be drawn with probability po and 
the probability of a coupon from collection (triple state) 
of size n is ^ i = 1,2,..., 24. 

Thus, we treat all 24 error states as one undesired 
coupon of probability p 0 = 24 x , where p is the prob¬ 
ability that the measurement setup is perfectly aligned. 
Let T n be a random variable describing a number of col¬ 
lected coupons until elements of every type are found for 
the first time. Then for n coupons with almost-uniform 
probability distribution and one extra element the aver¬ 
age value of T n is given by 


E (T„) 


n 


1 - Po 


E 


i 

i 


(15) 


We recalled this result from m , where far more general 
case was considered. For perfect measurements po = 0 
we get T n ft 90.6. 




